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ABSTRACT 

There are many Slowly Pulsating B (SPB) stars and 7 Dor stars in the Kepler Mission data set. 
The light curves of these pulsating stars have been classified phenomenologically into stars 
with symmetric light curves and with asymmetric light curves. In the same effective temper¬ 
ature ranges as the 7 Dor and SPB stars, there are variable stars with downward light curves 
that have been conjectured to be caused by spots. Among these phenomenological classes of 
stars, some show ‘frequency groups’ in their amplitude spectra that have not previously been 
understood. While it has been recognised that nonlinear pulsation gives rise to combination 
frequencies in a Fourier description of the light curves of these stars, such combination fre¬ 
quencies have been considered to be a only a minor constituent of the amplitude spectra. In 
this paper we unify the Fourier description of the light curves of these groups of stars, show¬ 
ing that many of them can be understood in terms of only a few base frequencies, which we 
attribute to g mode pulsations, and combination frequencies, where sometimes a very large 
number of combination frequencies dominate the amplitude spectra. The frequency groups 
seen in these stars are thus tremendously simplified. We show observationally that the com¬ 
bination frequencies can have amplitudes greater than the base frequency amplitudes, and we 
show theoretically how this arises. Thus for some 7 Dor and SPB stars combination frequen¬ 
cies can have the highest observed amplitudes. Among the B stars are pulsating Be stars that 
show emission lines in their spectra from occasional ejection of material into a circumstellar 
disk. Our analysis gives strong support to the understanding of these pulsating Be stars as 
rapidly rotating SPB stars, explained entirely by g mode pulsations. 

Key words: asteroseismology - stars: interiors - stars: oscillations - stars: variables - 
stars: individual (KIC 5450881; KIC7468196; KIC 8113425; KIC 10118750; KIC 10799291; 
Kiel 1971405) 


1 INTRODUCTION 

It is well known that Fourier analysis of non-sinusoidal light curves 
gives rise to harmonics and combination frequencies that describe 
the light curve shape in terms of sinusoids. High amplitude pul¬ 
sations are nonlinear, giving rise to significant amplitudes at the 
harmonics of the base frequencies. Multi-mode nonlinear pulsation 
results in interaction among the base frequencies and their harmon¬ 
ics that give rise to sum and difference combination frequencies 
of the form nfi ± mfj. Studying the relationships among the am¬ 
plitudes and phases of the Fourier components has been standard 
practice for RRLyr stars and Cepheids since the pioneering work 


of Simon & Lee (1981), and the study of combination frequencies 
and their astrophysical implications is well established for white 
dwarf stars (see, e.g. Wu 2001 and Montgomery 2005). Combina¬ 
tion frequencies dominate the amplitude spectra of some 5 Set stars, 
for example KIC 11754974 (Murphy et al. 2013) and KIC 8054146 
(Breger et al. 2012; Breger & Montgomery 2014), where the astro- 
physical implications and uses of the combination frequencies are 
more uncertain than for white dwarf stars. 

Among B, A and F main-sequence stars there are two classes 
of g mode pulsators: the 7 Dor stars with temperatures in the range 
of early to mid-F stars, and the Slowly Pulsating B (SPB) stars with 
temperatures in the range of the late B stars. The light curves of 
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these stars have been widely discussed phenomenologically, partic¬ 
ularly in the era of the photometric space missions MOST, CoRoT 
and Kepler. Debosscher et al. (2011) performed an automated vari¬ 
ability analysis on about 150 000 light curves from the Kepler Quar¬ 
ter 1 data, finding many 7 Dor stars. Balona et al. (2011a) visually 
scanned about 10000 stars in the Kepler data in the temperature 
range of the 7 Dor stars and the coolest SPB stars and classified 
the light curves phenomenologically as symmetric or asymmetric, 
where the asymmetric light curves show larger range at maximum 
amplitude than at minimum amplitude. Tkachenko et al. (2013) de¬ 
termined atmospheric parameters from high resolution spectra for 
69 stars in the Kepler data set that have 7 Dor g-mode pulsations. 
Balona et al. (2011b) similarly provided visual descriptions of 48 
B stars in the Kepler data. Following the lead of these papers, Mc¬ 
Namara, Jackiewicz & McKeever (2012) classified the light curves 
of 252 B stars in the Kepler data, describing many of the stars as 
‘Fg’, meaning that they show frequency groups in their amplitude 
spectra. Most recently Bradley et al. (2015) searched among 2768 
Kepler stars for 7 Dor stars, 5 Set stars and so-called ‘hybrid’ stars 
that show both p-mode and g-mode pulsations. They adopted the 
notation of Balona et al. (2011a) to describe the light curves as 
‘symmetric’ and ‘asymmetric’. 

Many of these papers used limited data sets from Kepler, un¬ 
intentionally resulting in significant confusion in the description of 
the light curves. In addition to pulsation, stars may show light vari¬ 
ability caused by orbital or rotational variations. Those are typically 
non-sinusoidal, hence also give rise to harmonics of the base fre¬ 
quencies. They do not, however, generate combination frequencies. 
It is thus possible to distinguish pulsation from rotational or orbital 
variability when combination frequencies are present. 

From the unprecedented time span of 4 yr of the full Kepler 
data set, it is now clear that g modes in 7 Dor stars and SPB 
stars can be so closely spaced in frequency that data sets span¬ 
ning less than 1 yr may not resolve the individual pulsation fre¬ 
quencies. Excellent examples of this are seen in the 7 Dor - 5 Set 
stars KIC 11145123 (Kurtz et al. 2014; Van Reeth et al. 2015) 
and KIC 9244992 (Saio et al. 2015), and in several other exam¬ 
ples given by Bedding et al. (2014) and Van Reeth et al. (2015), 
where the frequency spacings of long series of consecutive radial 
overtone g modes with rotational multiplets require up to half a year 
of data for full resolution. Thus previous descriptions of the char¬ 
acter of 7 Dor and SPB light curves based on relatively short data 
sets should be viewed with caution. The visual descriptions of light 
curves where pulsation modes are not resolved and combination 
frequencies are not recognised have led to erroneous conclusions. 

The presence of combination frequencies in 7 Dor and SPB 
stars has been recognised by many. Degroote et al. (2009) devel¬ 
oped an automated combination frequency search for CoRoT SPB 
stars. Papics (2012) discussed in general the search for combination 
frequencies for B, A and F main-sequence stars and problems asso¬ 
ciated with their identification, while Balona (2012) gave a detailed 
discussion in the case of 5 Set stars, particularly in comparison with 
white dwarf stars with the purpose of eventually using the informa¬ 
tion in the combination frequencies for astrophysical inference. 

Nevertheless, as Balona (2012) pointed out, the combination 
frequencies have usually been considered a nuisance in the search 
for pulsation mode frequencies for asteroseismology. He also states 
that ‘combination frequencies are of much lower amplitude than 
parent mode .. . [frequencies] . ..’. While this may be a widely held 
view, it is not necessarily true: Combination frequencies can have 
observed amplitudes greater than those of the base frequencies, as 
we show theoretically in Section 2. 


Another misconception is that a light curve that has larger 
variation at minimum than maximum cannot be purely pulsational. 
This idea has led to interpretations of some SPB and pulsating Be 
star light curves as being caused by spots, either completely or par¬ 
tially. The idea has crossed over into the visual description of the 
light curves of 7 Dor stars such that the papers using the termi¬ 
nology ‘asymmetric’ for the non-sinusoidal pulsators only include 
those stars that show more variation at maximum light than mini¬ 
mum light, even though there are 7 Dor stars that do the opposite, 
as is common for SPB stars. We show examples in sections below. 

We use our own description of the light curves at some ex¬ 
pense of proliferating nomenclature. We describe stars that have 
non-sinusoidal light curves with larger range at maximum light 
than minimum light as having ‘upward’ light curves, and those that 
do the opposite as having ‘downward’ light curves. Stars previously 
classified as having symmetric light curves are part of a continuum 
between these extremes. Below we show examples of the various 
shapes of the light curves and their simple explanation in terms of 
non-sinusoidal pulsation in only a few pulsation modes with com¬ 
bination frequencies. We find that many 7 Dor, SPB and pulsating 
Be star light curves are far simpler than has previously been under¬ 
stood, and we make a strong case that the only physics needed to 
understand all of these stars is nonlinear pulsation theory. There is 
no need of, and no evidence for, spots. 

The reduction that we demonstrate in the apparent complex¬ 
ity of the amplitude spectra of the stars showing frequency groups 
is stunning. Instead of hundreds of frequencies being extracted 
for analysis, many of these stars have but a few pulsation mode 
frequencies with a plethora of combination frequencies, some of 
which can have amplitudes greater than the base frequencies. As in 
the cases of the p-mode pulsation in the 5 Set stars KIC 11754974 
and KIC 8054146 mentioned above, the amplitude spectra of the 
g-mode pulsators on the main-sequence can be dominated by com¬ 
bination frequencies. These must be fully modelled to get to the 
pulsation mode frequencies that are the fundamental data of aster¬ 
oseismology, and they have the potential to provide new astrophys¬ 
ical information for main-sequence stars, as they do for pulsating 
white dwarf stars. 

In this paper we explain that what previously appeared to be 
complex variability with dozens or hundreds of frequencies is a re¬ 
sult of only a few base frequencies and their combination frequen¬ 
cies. This is an important observational result that greatly simplifies 
our understanding of the light curves of 7 Dor, SPB and pulsat¬ 
ing Be stars. That such strong nonlinear interaction exists indicates 
high amplitudes for the base modes in the stellar cores. It is a goal to 
gain asteroseismic inference from these modes by modelling them. 


2 THEORETICAL CORROBORATION 
2.1 Weak nonlinear system 

To draw some basic characteristics of nonlinear pulsation of stars, 
let us consider a case of weakly nonlinear pulsation, in which the 
eigenfrequencies are still close to those obtained by linear calcula¬ 
tion. We consider the unperturbed static equilibrium state of a star 
and superimpose on it perturbations. To make the problem simple, 
we assume that in the equilibrium state the star is spherically sym¬ 
metric. 

We define the displacement vector, 
i{ro,t) := r - ro, ( 1 ) 

where r denotes the Lagrangian position variable of a given fluid 
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element which is at r = ro in the equilibrium state. The equation 
of oscillations, which is expressed with a single variable is then 
divided into the linear operator C{$,) and the nonlinear operators 

(fe = 

^+£{$.)+ (^, (^, ^) + • • • = 0 , ( 2 ) 

where denotes the operator of the fc-th order of Retaining 
only the first-order terms, we obtain the linearised equation 

^+^(0=0. (3) 

Since £ does not include any operator with respect to time, the 
solution to equation (3), defined as is separated into a 

spatial function and a temporal function. The latter is expressed by 
exp iujt, with frequency w. Equation (3) turns into an eigenvalue 
problem with a set of suitable boundary conditions. The eigenfunc¬ 
tions form an orthogonal complete set, hence the linear adiabatic 
oscillations of a star can be expressed as 

= '^ah'E,^^\ro)expi{iUktifik), (4) 

fc 

where tUfc and H(ro) denotes the eigenfrequency and the eigen¬ 
function of the mode index fc, respectively, and Ofe and ipk are the 
amplitude and the phase of the mode at f = 0. Here, the mode in¬ 
dex fc consists of the spherical degree Z, the azimuthal order m, and 
the radial order n, and the eigenfunction is written with 

the spherical coordinates (r, 0,(j)) as 


= En,i(r)Yr Sr + 



1 ay,™ 

sin 6 dcj) ^ 


(5) 


where (j>) denotes the spherical harmonics with the spheri¬ 

cal degree I and the azimuthal order m, and E„,i{r) and Hn,i{r) 
are the radial eigenfunctions, with respect to r, for the displace¬ 
ment in the radial direction and for that in the horizontal direction, 
respectively, with the radial order n and the spherical degree 1. The 
characteristics of the linear adiabatic pulsations of stars have al¬ 
ready been investigated in detail, as in the textbooks Unno et al. 
(1989) and Aerts, Christensen-Dalsgaard & Kurtz (2010). 

Including the second-order terms of equation (2), we obtain 
the following equation for t): 


^ 2 ^( 2 ) 

dd 




( 6 ) 


Since has already been independently solved, the 

above equation (6) is regarded as an inhomogeneous equation for 
with a source term which originates 

from the squared terms of the first-order free oscillations. In other 
words, equation (6) is regarded as an equation for a forced oscilla¬ 
tion induced by the nonlinear term The particu¬ 

lar solution to this inhomogeneous equation gives the correction to 

In a similar way, the higher-order solutions are considered as 
forced oscillations successively induced by the nonlinear terms of 
the lower-order solutions. 


2.2 Why do the comhination frequencies appear? 

The operator consists of cross terms of (ro, Z), and it is 
bilinear. Hence, the nonlinear term is separated into a spatial part 


and a temporal function. Substitution of the form of given by 
equation (4) into equation (6) leads to 


dd 




fc,fc' 








The cross terms in the nonlinear operator induce the com¬ 
bination frequencies. As a consequence, the particular solution to 
(ro, t) also has combination frequencies -f 
It should be noted that the associated general solution of the 
inhomogeneous differential equation (7) is of the form of equation 
(4), and it is already given as a first-order solution. Hence we only 
have to consider the particular solution. 

A special case is the cross term of u>k with itself. That induces 
the second harmonic 2uJk- Similarly, the nonlinear operator A/”*-^^ 
produces the third harmonic 3cJfc through the cross term between 
ujk and 2cJfc or the triple term of ujk- This is the process producing 
a non-sinusoidal light curve from a single mode. 


2.3 Why do some combination frequencies have amplitudes 
greater than their base frequencies? 

The second-order perturbation is of the order of the square of the 
linear perturbation. It should be noted here, however, that this state¬ 
ment concerns the intrinsic amplitudes. The visibility, which is 
highly dependent on the surface pattern of the oscillations, must 
be taken into account to evaluate the actual observed amplitudes. 

The nonlinear operator Af induces cross terms of spherical 
harmonics, and they are described in terms of a series of spheri¬ 
cal harmonics with azimuthal order that is equal to the sum of the 
parent spherical harmonics; 

Yr{eA)yf{e,<j>) 

= " (Z, -m, z', m')y, [e, 0), (8) 

l" 

where —m, I', m!) is defined by 

:= J dnYr{e,<j>rYf{e,<i>)Yr"\e,<j>),(9) 

and I” is in the range of [ | Z — Z' |, Z -|- Z' ], except the range of 
[0,1 m -b m' I ]. This means that even if the first-order perturba¬ 
tions associated with high degree have low observed amplitudes 
and are difficult to detect, their products may induce low degree 
components, e.g., Z = 0, that are detectable. 

So, it is not necessarily true that combination frequencies of 
higher-order perturbations have smaller observed amplitudes than 
the base frequencies. Some combination frequencies can have ob¬ 
served amplitudes greater than those of their base frequencies. 


2.4 Why do some stars show downward light curves? 

Superposition of two oscillations with nearly equal frequencies 
leads to a beat phenomenon. With an increase in the frequency 
difference, the wave pattern gradually changes. When the second 
harmonic is imposed on the base frequency, the oscillation pattern 
significantly deviates from symmetry with respect to the zero level. 
If the phase difference of these two frequencies, iph — Pk', is nearly 
zero, the wave pattern shows an ‘upward’ shape (which has previ¬ 
ously been called ‘asymmetric’), whereas when the phase differ- 
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time 

Figure 1. This simulated light curve demonstrates a combination of fre¬ 
quencies differing by a factor of two with a phase difference of 7t induces a 
‘downward’ light curve. 


Table 1. Basic data for the stars presented in this paper. Teff and log g are 
from Huber et al. (2014). The last column gives the section of the paper 
where the star is discussed. 


KIC 

Teff 

K 

logs 

cgs units 

Kp 

mag 

type 

Section 

5450881 

10500 ± 250 

4.1 ± 0.2 

12.45 

SPB 

6.3 

7468196 

6850 ± 200 

4.1 ± 0.2 

13.77 

7 Dor 

5 

8113425 

6900 ± 200 

4.3 ± 0.1 

13.86 

7 Dor 

4 

10118750 

11400 ± 400 

3.7 ± 0.3 

13.90 

SPB 

6.1 

10799291 

11100 ±400 

4.4 ± 0.1 

14.98 

SPB 

5 

11971405 

11600 ±400 

3.7 ± 0.4 

9.32 

SPB/Be 

7.1 


ence is close to n, the wave pattern shows a ‘downward’ shape as 
shown in Fig. 1. 

It has been widely considered that nonlinear pulsation induces 
only upward light curves and that pulsation cannot be the sole phys¬ 
ical cause of downward light curves. However, as demonstrated 
here, pulsation induces downward light curves in some cases. The 
base frequencies are driven by the K-mechanism or by convective 
blocking, while the harmonic frequencies are damped by heat loss. 

Differences in thermal properties may cause phase differences 
between these two extreme groups to differ by tt, with intermediate 
phases giving rise to less extreme distortion of the light curves. We 
thus propose that the range of pulsational light curves in B, A and F 
main-sequence stars, from upward through symmetric to downward 
shapes, is a consequence of the phases of the nonlinear harmonic 
and combination frequencies, and that those phases are determined 
by the balance between driving and damping in each individual star. 


3 DATA AND ANALYSIS METHODS 

We have visually examined Kepler light curves and amplitude spec¬ 
tra for thousands of B, A and F main-sequence stars and selected 
examples to illustrate our results. The data used for the analysis in 
this paper are the Kepler quarters 0 to 17 (QO - Q17) long cadence 
(LC) data. The Kepler ‘quarters’ were of variable time span that de¬ 
pended on operational constraints. Most quarters are close to one 
fourth of a Kepler orbital period of 372.455 d, which was the time 
scale on which the satellite was rolled to keep the solar panels fully 


illuminated. QO, Q1 and Q17 were short ‘quarters’. We do not use 
any stars in this paper that fell on the failed module 3, so there are 
no large gaps in our data sets. 

We used the multi-scale, maximum a posteriori (msMAP) 
pipeline data converted to magnitudes; information on the reduc¬ 
tion pipeline can be found in the data release notes^ 21. To op¬ 
timise the search for exoplanet transit signals, the msMAP data 
pipeline removes some astrophysical signals with frequencies less 
than 0.1 d“^ (or periods greater than 10 d). Some of the combina¬ 
tion frequencies that we find in this paper are at frequencies less 
than this 0.1 d“^ limit, and we show that these frequencies are un¬ 
perturbed by the pipeline reductions. A useful general conclusion is 
that while the pipeline may reduce astrophysical amplitude, it does 
not perturb the frequencies. 

For all stars we visually inspected the light curve and removed 
a few obvious outliers. The time span for QO-17 data is 1470.5 d, 
and for Q1-Q17 data is 1459.5 d. Typically, about 64 840 data 
points comprised each data set. Table 1 gives basic data for the stars 
analysed in this paper. 


3.1 Frequency analysis 

We first produced a catalogue of all stars in the Kepler data with 
effective temperatures in the Kepler Input Catalogue (KIC) above 
6400 K. Our catalogue comprises light curves and amplitude spec¬ 
tra for each quarter of Kepler data for each star. We have visually 
studied each of these plots. For stars that we studied in more detail, 
we first examined the entire data sets, Q0-Q17, using the interac¬ 
tive light curve and amplitude spectrum tools in the programme 
PERIOD04 (Lenz & Breger 2004). We then used a Discrete Fourier 
Transform (Kurtz 1985) and our own least-squares and nonlinear 
least-squares fitting programmes to find the frequencies, ampli¬ 
tudes and phases to describe the light curves. After normalising 
the entire data set to zero in the mean, we fitted a cosine function. 
Am = Acos(27r/(f — to) + 4>), to the data in magnitudes, thus 
defining our convention for the phases in this paper. Our routines 
and PERIOD04 are in agreement. 

Our procedure was to identify ‘base frequencies’ in the ampli¬ 
tude spectrum from which to generate the combination frequen¬ 
cies. Those were then optimised by fitting them simultaneously 
by non-linear least-squares to the data. For reasons of space and 
presentation we do not tabulate the individual frequency uncer¬ 
tainties. Those depend on the signal-to-noise ratio (Montgomery 
& O’Donoghue 1999) and are in general of the order of 10~^ — 
10“®d“^. This is significantly less than the resolution of the 
Fourier peaks of 77 ~ 1/AT = 0.0007d“^, where T = 1460d 
is the time span of the data, but the frequencies are much better 
determined than the resolution of the data set, so long as there are 
no unresolved frequencies. We discuss both cases in the sections 
below. 

Importantly, following the determination of the base frequen¬ 
cies, we did not extract large numbers of peaks in the amplitude 
spectra, then test whether peaks were combination of the base fre¬ 
quencies. It is so clear from first inspection that the amplitude spec¬ 
tra of many 7 Dor and SPB stars are dominated by combination 
frequencies that we selected just a few base frequencies and then 
calculated the frequencies of the combination terms and fitted that 
calculated set of frequencies by least-squares to the data. The suc- 


^ https://archive.stsci.edu/kepler/data_release.html 
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Figure 2. A section of the light curve of KIC 8113425 spanning 100 d show¬ 
ing the strongly nonlinear ‘upward’ light variations. The time is relative to 
B JD 245 5000. Almost all of the variation in the light curve is explained 
by the nonlinear interaction of only four g modes. 


cess of that procedure in accounting for most of the variance in the 
data sets is apparent in the sections below and justifies the method. 

Of course, because of the nature of combination frequencies, 
other sets of base frequencies chosen from the combination fre¬ 
quencies can produce the same set of frequencies to be fitted to 
the data. Hence the identification of base frequencies with an as- 
trophysical cause is open to interpretation. We hypothesise in this 
paper that all of the chosen base frequencies arise from g mode pul¬ 
sations. Other interpretations, should they arise, do not change the 
mathematical fit of the chosen cosinusoids to the data. 


4 THE 7 DOR STAR KIC 8113425 

Frequency groups were noted in some Kepler stars by McNa¬ 
mara, Jackiewicz & McKeever (2012) following the terminology 
of Balona et al. (2011a), which we adopt. Our own examination of 
the light curves and amplitude spectra of thousands of Kepler B, A 
and F stars shows that many Kepler 7 Dor and SPB stars have such 
frequency groups with a wide variety of characteristics, as has been 
noted by others. For example, some illustration of this can also be 
seen in Figure 2 of Tkachenko et al. (2013) who studied 69 7 Dor 
stars in the Kepler data set. 

These frequency groups may arise from different causes, but 
many of the stars showing them have only a few nonlinear pulsa¬ 
tion modes with amplitude spectra that are dominated by combi¬ 
nation frequency peaks. We illustrate this with an extreme case: 
KIC 8113425 is a 7 Dor star with an upward light curve, which is 
shown in Fig. 2 with a 100-d section of the light curve. The strongly 
nonlinear pulsations can be seen clearly. This star previously was 
listed as an ‘asymmetric’ light curve star that shows ‘evidence of 
migrating starspots’ according to Balona et al. (2011a). We argue 
here that the star has no spots; that suggestion came from a sim¬ 
ple visual inspection of a short section of the light curve with no 
frequency analysis. 

KIC 8113425 was observed by Kepler in LC in all quarters 
from Q1 to Q17. No other published observations or data are avail¬ 
able. Fig. 3 shows an amplitude spectrum of the Ql-17 Kepler data 
out to a frequency of 2 d~^ where it can be seen that there are many 
frequency groups, which we label fgO, fgl, etc. In a close exami¬ 
nation of the amplitude spectrum, frequency groups up to fgl 1 can 
be seen, showing the presence of very high order combination fre¬ 
quencies. 

We chose KIC 8113425 to discuss the understanding of the 
frequency groups because it shows so many groups which are well- 
separated in frequency. It appears to be complex and to be an ex¬ 
treme example of the frequency group phenomenon, yet the appar¬ 
ent complexity lies almost entirely in the combination frequencies 
- we have fitted most of the variance with only four base frequen¬ 
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Figure 3. Top: An amplitude spectrum of the Kepler Ql-17 data for 
KIC 8113425 out to 2d~^. There are no p mode pulsations at higher fre¬ 
quencies up to the Nyquist frequency. The frequency groups are labelled up 
to fg4, but can be seen to extend further; at least 11 groups can be seen in the 
amplitude spectrum of this star. The middle panel shows on the same scale 
an amplitude spectrum of the residuals after pre-whitening the four pulsa¬ 
tions frequencies and 39 combination frequencies with terms up to order 2v 
given in Table 2, and with amplitudes greater than 1 mmag; the inclusion of 
combination frequencies with amplitudes less than 1 mmag removes more 
peaks. Higher order combination frequencies explain the higher frequency 
groups. The reduction in variance is stunning. The four base frequencies 
from frequency group 1 (fgl) are shown and labelled in the bottom panel. 
We propose that those represent g mode pulsations. 


cies. Table 2 lists those four base frequencies and 39 combination 
frequencies with terms up to order 2v that were fitted to produce 
the amplitude spectrum of the residuals shown in the middle panel 
of Fig 3. 

We limit our discussion here to combination frequencies with 
terms to order 2v and amplitudes greater than 1 mmag so that we 
have a manageable number of frequencies from which to discuss 
the details. A fit of combination frequencies with terms up to order 
bu with no lower limit on amplitude yields over 500 combination 
frequencies and reduces the variance in the amplitude spectrum of 
the residuals further. Obviously, to explain the highest frequency 
groups would need even higher order combination frequencies. As 
Papics (2012) warned, with higher order fits of combination fre¬ 
quencies, the number of fitted peaks can become large with respect 
to the number of independent Fourier peaks in the amplitude spec¬ 
trum. Thus in future detailed studies of these stars with frequency 
groups, care will need to be taken with high-order fits to avoid 
combination frequencies that occur only by coincidence. With the 
large amplitudes and low number of combination frequencies we 
use here, such chance coincidences are not a problem. 
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Table 2. A least-squaies fit of the four pulsation mode frequencies of 
KIC 8113425 and their combination frequencies with terms up to order 2v 
having amplitudes greater than 1 mmag. There are 43 identified frequen¬ 
cies, including the four base frequencies. The zero point of the time scale is 
BJD245 5694.25. 


labels 

frequency 

d-i 

amplitude 

mmag 

±0.03 

phase 

radians 

fgo 

— ui + 21/3 - 1/4 

0.003055 

1.87 

-1.612 ± 0.014 

Ui — U 2 — Us + i /4 

0.008108 

4.18 

-1.757 ± 0.007 

-U2 + i/3 

0.011163 

1.20 

1.088 ±0.025 

-i /3 + U 4 , 

0.028150 

3.71 

-1.073 ±0.008 

-i/r + Us 

0.031206 

2.42 

1.997 ± 0.012 

-i/l + Ui 

0.059356 

4.32 

0.874 ±0.007 

- 2 ui + US + U 4 

0.090562 

1.37 

-0.033 ±0.021 

fgl 

to 

1 

0.398853 

3.66 

1.299 ±0.008 

to 

1 

0.410016 

1.02 

2.107 ± 0.030 

U 2 + Us - Ui 

0.421951 

1.04 

2.913 ±0.029 


0.430058 

14.55 

0.742 ± 0.002 

^2 

0.450101 

10.09 

-0.449 ± 0.003 

UI - US + Ui 

0.458209 

1.95 

-2.793 ±0.015 

US 

0.461264 

7.75 

-1.995 ± 0.004 

— Ul + U 2 +US 

0.481307 

1.73 

1.891 ± 0.017 

1^4 

0.489414 

13.14 

-3.136 ±0.002 

— U 2 + US +Ui 

0.500577 

1.57 

-0.331 ± 0.019 

-Ul + U 2 +Ui 

0.509457 

1.47 

-2.478 ±0.021 

— Ul + US +Ui 

0.520620 

1.51 

2.294 ± 0.020 

fg 2 

2ui 

0.860117 

2.42 

-1.921 ± 0.012 

Ul + U 2 

0.880159 

1.21 

-2.729 ±0.025 

2ui - i/3 + i/4 

0.888267 

1.91 

1.134 ±0.016 

Ul + i /3 

0.891322 

4.79 

1.129 ±0.006 

Ul + Ui 

0.919473 

7.39 

0.647 ± 0.004 

2us 

0.922528 

2.16 

-2.381 ± 0.014 

Ul — Us + 2ui 

0.947623 

1.29 

-2.936 ±0.023 

Us + Ui 

0.950678 

6.48 

-2.749 ± 0.005 

2ui 

0.978828 

4.10 

-3.036 ± 0.007 

-Ul + 2us + Ui 

0.981884 

1.01 

2.974 ± 0.029 

-Ul + i /3 + 2ui 

1.010034 

1.63 

2.264 ±0.018 

fg3 

2,ui + Us 

1.321381 

1.50 

-1.083 ±0.020 

2ui + U 4 

1.349531 

1.64 

-1.520 ±0.018 

Ul + 2 us 

1.352586 

1.02 

2.184 ±0.029 

Ul +U 2 + Ui 

1.369574 

1.02 

-2.624 ±0.029 

2ui — i/'s + 2v4 

1.377681 

1.15 

1.409 ±0.026 

Ul + us + Ui 

1.380736 

3.33 

1.522 ± 0.009 

U2 + Us + Ui 

1.400779 

1.45 

0.119 ±0.021 

Vi + 2 i /4 

1.408887 

1.88 

0.814 ±0.016 

2u^ + V4 

1.411942 

1.04 

-1.284 ±0.029 

US + 2 ui 

1.440092 

1.07 

-1.964 ±0.028 

fg4 

2i/l + i /3 + i /4 

1.810795 

1.87 

-0.920 ±0.016 

Ul + 2us + Ui 

1.842000 

1.61 

2.097 ± 0.018 

Ul +US + 2ui 

1.870151 

1.94 

1.606 ±0.015 



Time (BJD - 2455835 d) 



Figure 4. Top: Kepler observations of KIC 81 13425. Bottom: an artificial 
light curve (black) constructed with a single frequency [red = cos(27ra)] 
and its harmonic [blue = 0.5cos(2 X 27ra;)]. The relative phase of the 
harmonic is zero, so it reinforces the maxima and suppresses the minima, 
leading to an ‘upward’ light curve. Similarities with the observations in the 
top panel are already evident after including only one combination term. 

4.1 The phases 

The asymmetry of the light curve is described through the combi¬ 
nation frequencies, but the contribution of each frequency to that 
description differs. Three factors are important in determining the 
asymmetry described by each combination. We explain them here 
for KIC 8113425, where the upward light curve shows high max¬ 
ima and shallow minima, but the description could also be applied 
to a downward light curve with low maxima and deep minima, or 
any intermediate shape light curve. The discussion can thus be gen¬ 
eralised from this specific example. 

The first important factor is the combination frequency. In or¬ 
der to describe sharp and high maxima, the combination frequency 
must be higher than the base frequency. The second factor is the 
phase. To describe high maxima in the light curve, the combina¬ 
tions need a relative phase close to zero. The relative phase^, (pr, of 
a combination frequency is defined as 

(pT — (pohs (5calc — (pohs “t“ ) (10) 

where c^oba is the observed phase of the combination frequency, 
and (pcaic is a phase calculated from the base frequencies, which in 
this case is for the combination v = nvi -|- mvj. A relative phase 
of zero means the maximum of the combination frequency coin¬ 
cides with the maxima of the base frequencies, and this maximum 
is then reinforced. An example is shown in Fig. 4. Finally, the third 
factor is amplitude. Combinations with low amplitudes have little 
relevance to the shape of the light curve. 

The contribution each combination frequency makes to the de¬ 
scription of the asymmetry can be shown in a phasor (a contraction 
of ‘phase vector’) diagram. Fig. 5. Here, the amplitude and relative 
phase of each combination are shown on a polar plot. The conven¬ 
tions chosen dictate the orientation of the phasor diagram; here we 

^ We calculated the relative phases with a cosine function applied to the lu¬ 
minosity variations so that upward light curves have combination frequency 
phases in the positive x direction in the phasor plots, whereas, for purposes 
of pre-whitening in the frequency analysis where logarithms are preferred, 
the phases in the tables are for magnitude variations. There is thus a Tirad 
shift in each phase in the tables that was used to calculate the phases in the 
phasor plots. 
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Figure 5. A phasor plot for KIC 8113425 showing the relationship of the 
phases of the combination frequencies to those of the base frequencies 
(black stars), which by definition have relative phase zero (i.e., they lie 
along the positive Cartesian x-axis, which represents (pr = 0). Variables 
shown are amplitude and relative phase of the combination frequencies, 
using cosines to describe the luminosity variation in units of parts-per- 
thousand. Hence, points that ai‘e right of centre correspond to combination 
frequencies describing ‘upward’ asymmetry. The groups described with dif¬ 
ferent colours and symbols correspond to the frequency groups of Table 2. 


have chosen cosines to fit the luminosity variation, so that points 
that lie to the right of centre correspond to combination frequen¬ 
cies describing ‘upward’ asymmetry. Similarly, points to the left of 
centre belong to combinations describing ‘downward’ asymmetry. 
Since KIC 8113425 has a strong upward asymmetry, we see points 
of high amplitude in the right part of the diagram. 
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Figure 6. Top panel: A section of the light curve of the 7 Dor star 
KIC 7468196 spanning 40 d showing the upward light variations. Bottom 
panel: A section of the light curve of the SPB star KIC 10799291 spanning 
the same 40 d showing the downward light variations. The time is relative 
to BJD 245 5000. Note that while the light curves have a different appear¬ 
ance, they are both explained by a few nonlinear g-mode pulsations with 
combination frequencies. 


5 A 7 DOR STAR AND AN SPB STAR: UPWARD AND 
DOWNWARD LIGHT CURVES COMPARED 

As for 7 Dor stars such as KIC 8113425 shown in Section 4, many 
SPB stars also show frequency groups that have previously lacked 
an explanation (e.g., McNamara, Jackiewicz & McKeever 2012). 
In this section we compare the light curves and amplitude spectra 
for a 7 Dor star and an SPB star, showing a strong similarity that 
demonstrates that the SPB star light curve is also fully explained 
by a few base frequencies with combination frequencies, just as is 
the 7 Dor star light curve. The 7 Dor star, KIC 7468196, has an 
upward light curve, and the SPB star, KIC 10799291, has a down¬ 
ward light curve. We show that this difference in light curve shape 
arises only from the phases of the combination frequencies, and 
we conclude that this shows that the SPB star, with its downward 
light curve, is fully explained by a few g-mode pulsation frequen¬ 
cies and their combination frequencies. Fig. 6 compares typical 40- 
d sections from the light curves of the two stars. 

Fig. 7 compares the amplitude spectra for the two stars, show¬ 
ing remarkable similarities. Both amplitude spectra are the result 
of a few base frequencies plus combination frequencies. For the 
7 Dor star, KIC 7468196, a frequency septuple! can be seen in the 
frequency range labelled as fg2 in the figure. Another striking fre¬ 
quency septuple! is seen for the SPB star, KIC 10799291, also in 
the frequency range fg2. In both cases the central frequency of the 
fg2 septuple! is the sum of two of the base frequencies, and the 
other components are further combination frequencies. It is partic¬ 
ularly notable that for KIC 10799291 the amplitude of the high¬ 
est peak in fg2, 1/2 -F 1^7 in our notation, has a higher amplitude 
than either of its two base frequencies. For pulsating stars of the 
main-sequence it has not previously been recognised that the com¬ 
bination frequency amplitudes can be this large, hence they have 
usually not been recognised. We derived theoretically in Section 2 
that this is possible, and we demonstrate here observationally that 
it happens. 

For the 7 Dor star KIC 7468196 we fitted only three base fre¬ 
quencies and their combination frequencies with terms up to order 
3zz, with the results shown in Table 3. As can be seen in the mid¬ 
dle left panel of Fig. 7, the reduction in the variance is striking. For 


Table 3. A least-squares fit of the three pulsation mode frequencies of 
the 7 Dor star KIC 7468196 and the combination frequencies with terms 
up to order for base frequencies ui, U 2 and vz and with amplitudes 
restricted to greater than 0.3mmag. There are 27 identified frequencies, 
including the three base frequencies. The zero point of the time scale is 
BJD 245 5694.25 


labels 

frequency 

d-i 

amplitude 

mmag 

±0.004 

phase 

radians 

fgO 

— Ul + V2 

0.1012839 

4.556 

0.8699 ± 0.0009 

—2i^l -t- 21^2 

0.2025678 

1.890 

-0.4663 ± 0.0021 

— 1^2 + 173 

0.2337767 

1.190 

-0.9808 ± 0.0034 

-3i/i + 3 u2 

0.3038517 

0.678 

3.1209 ±0.0059 

-u\ + 173 

0.3350606 

1.993 

1.3341 ± 0.0020 

— 2 ui + U 2 + 173 

0.4363445 

0.606 

0.4556 ± 0.0066 

fgl 

3ui — 2 u 2 

1.1632242 

0.688 

1.1930 ±0.0058 

2 i 72 — U 3 

1.2332992 

1.520 

2.7069 ± 0.0026 

2 u\ — V 2 

1.2645081 

5.628 

-1.5188 ±0.0007 

+ 3u2 — 173 

1.3345831 

1.372 

1.8239 ±0.0029 

171 

1.3657920 

23.879 

-2.1298 ±0.0002 

3i^i — 3 v 2 + 173 

1.3970009 

0.333 

2.1111 ± 0.0120 

172 

1.4670759 

20.357 

0.6546 ± 0.0002 

2 ui — 2 v 2 + 173 

1.4982848 

0.392 

-0.2886 ±0.0102 

— Ui+ 2 v 2 

1.5683598 

1.778 

-1.3550 ±0.0023 

— 2 vi + 3172 

1.6696437 

0.988 

-2.4776 ± 0.0041 

173 

1.7008526 

5.319 

1.8704 ± 0.0008 

-171 +172 + 173 

1.8021365 

0.558 

-0.1811 ± 0.0072 

—2ui + 2 v 2 + 173 

1.9034204 

0.303 

-1.6220 ±0.0132 

fg2 

U 1 + 2 v 2 — 173 

2.5990912 

0.475 

-2.2965 ± 0.0084 

3ui — V 2 

2.6303001 

1.519 

-0.2981 ± 0.0026 

2 vi 

2.7315840 

2.066 

-0.8022 ±0.0019 

17l + 172 

2.8328679 

3.788 

1.9834 ± 0.0011 

2v2 

2.9341518 

2.220 

-1.4248 ±0.0018 

2ui - V 2 + V 3 

2.9653607 

0.323 

-2.6336 ±0.0124 

vi + 173 

3.0666446 

0.953 

-3.0223 ± 0.0042 

>72 + 173 

3.1679285 

0.686 

0.0280 ± 0.0058 


the SPB star KIC 10799291 there are more pulsation frequencies, 
and they are more widely spread in frequency, as can be seen in the 
right column of Fig. 7. Because of their relatively high amplitudes, 
we fitted eight pulsation frequencies, but only used four of those as 
base frequencies for calculating the combination frequencies with 
terms up to order 2v. The results are given in Table 4. 


5.1 The phases 

The contrasting shapes of the light-curve asymmetries in 
KIC 7468196 and KIC 10799291 lead to contrasting phasor plots, 
as shown in Fig. 8. The combination frequencies belonging to 
KIC 10799291 lie almost entirely on the left of the diagram, de¬ 
scribing strong, downward asymmetry (bottom panel). Conversely, 
the distribution of points in the upper panel is not strongly skewed 
to one side, indicating that the asymmetry in the light curve is not 
strong. Indeed, close inspection of the top panel of Fig. 6 shows that 
this is the case - the light curve has only mild upward asymmetry. 
The fact that we observe asymmetry at all, when the distribution is 
so even, is an effect of the frequencies of the combinations. Those 
that lie in fg2 have higher frequencies, and thus simultaneously 
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Figure 7. Left column: amplitude spectra for the 7 Dor star KIC7468196. Right column: amplitude spectra for the SPB star KIC 10799291. Top panels: 
Amplitude spectra of the Kepler Ql-17 data out to 4d“^ for both stars. There are no p-mode pulsations at higher frequencies up to the Nyquist frequency. 
Eye-catching features of these amplitude spectra are the frequency septuplets that make up the second frequency group, fg2, for both stars, but particularly 
for the SPB star on the right. These are consequences of combination frequencies of only three base frequencies. Middle panels: the amplitude spectra of the 
residuals after pre-whitening by the pulsation frequencies plus the combination frequencies of a subset of the base frequencies. These are given in Tables 3 
and 4. The scale is the same as in the top row, showing that almost all of the variance is explained for both stars. Bottom panels: Higher resolution looks at the 
fgl frequency ranges with pulsation mode frequencies labelled. In the case of KIC 10799291, only four of these were used as base frequencies to generate the 
combination frequencies, as shown in the tables. 


suppress the minima while reinforcing the maxima. The small up¬ 
ward asymmetry is thus explained by the small imbalance of the 
fg2 frequencies to the right of the diagram. 

5.2 Discussion of the comparison of KIC 7468196 and 
KIC10799291 

We argue from the similarities of the amplitude spectra of the 7 Dor 
star KIC 7468196 and the SPB star KIC 10799291, and from the 
fact that most of the variance in both light curves is explained hy a 
few base frequencies and their combination frequencies, that both 
stars are pulsating in g modes that fully account for their light vari¬ 
ations. While there are further details in the amplitude spectra at 
lower amplitudes that we have not extracted, we expect that these 
details will be accounted for by pulsation in other modes and the 
additional combination frequencies that these may generate. The 
difference in the appearance of the light curves, with the 7 Dor 
light curve being upward and the SPB star light curve being down¬ 
ward, is a consequence of the phases of the combination frequen¬ 
cies. There are 7 Dor stars with downward light curves, too. There 
is no need to conjecture spots to explain this, and spots do not gen¬ 
erate combination frequencies, hence are not a viable hypothesis. 
The light curves of the 7 Dor and SPB stars are fully explained by 
g mode pulsation. We will expand on this point in the final discus¬ 
sion. 

There has been some confusion about the classification of 
KIC I079929I that we clarify here. KIC I079929I has revised Ke¬ 
pler parameters Teff = 11100 ± 400 K and log g — 4.4 ±0.1 (cgs 


units) (Huber et al. 2014), as listed in Table 1. However, the origi¬ 
nal KIC parameters are T^s = 10950 K and logg = 6.1, which 
suggested a white dwarf star. McNamara, Jackiewicz & McKeever 
(2012) also classified this star as a white dwarf based on proper 
motion and the original KIC parameters. That classification propa¬ 
gated into SIMBAD, where KIC 10799291 is listed as a ‘pulsating 
white dwarf’. The revised KIC parameters of Huber et al. (2014) 
and the frequencies of the variation in the star show that it is an 
SPB star, not a white dwarf. 
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Table 4. A least-squares fit of the eight pulsation mode frequencies of 
the SPB star KIC 10799291 and the combination frequencies with terms 
up to order 2 ^, for base frequencies ui, U 2 , 1^3 and uj, and with ampli¬ 
tudes restricted to greater than 50 p-mag. There are 41 identified frequen¬ 
cies, including the eight base frequencies. The zero point of the time scale 
isBJD 245 5694.25 


labels 

frequency 

d-i 

amplitude 

mmag 

±0.004 

phase 

radians 

fgO 

2 ui + — 2 i /'7 

0.1357140 

0.047 

-2.292 ± 0.093 

— h'l + V 2 + 21/3 — 2 uj 

0.1462140 

0.096 

2.126 ± 0.045 

-U 3 + U 7 

0.1520070 

0.221 

1.719 ± 0.020 

— U 2 + 1^7 

0.1734908 

1.297 

-0.322 ± 0.003 

— VI + 2 v^ — V 7 

0.3197048 

0.039 

1.992 ± 0.112 

— 2 v 2 + 2^*7 

0.3469816 

0.106 

-1.748 ± 0.041 

— vi + 2 v 2 — V 2 

0.4287442 

0.053 

2.803 ± 0.082 

— vi + V 2 

0.4502280 

0.274 

1.069 ± 0.016 

2 vi — 2 iv 3 + V7 

0.5917350 

0.043 

-0.009 ± 0.101 

— VI + V 2 — 1^3 + I'T 

0.6022350 

0.263 

0.251 ± 0.017 

— Vl + VJ 

0.6237188 

0.426 

0.909 ± 0.010 

vi + V 2 — vj 

0.7379490 

0.049 

-0.582 ± 0.088 

fgi 


0.9114398 

0.565 

-1.668 ± 0.008 

— 2 vi + 2 i /'3 

0.9434236 

0.188 

-1.837 ± 0.023 

1^1 — 1^2 + 1^7 

1.0849306 

0.082 

-1.998 ± 0.053 

2 v 2 + 1^7 

1.1881770 

0.160 

-0.889 ± 0.027 

^2 

1.3616678 

2.474 

-1.487 ± 0.002 

^3 

1.3831516 

0.841 

-3.141 ± 0.005 

V4 

1.4214710 

0.286 

0.616 ± 0.015 


1.4754643 

0.671 

-0.629 ± 0.007 

2vi + 2v2 — 2 i-'7 

1.4758980 

0.022 

2.876 ± 0.222 

— Vl + 2v2 -\- 2 i/'3 — 2 i/'7 

1.5078818 

0.111 

0.317 ± 0.067 


1.5081214 

0.449 

1.902 ± 0.017 

t'7 

1.5351586 

1.927 

-1.626 ± 0.002 

^8 

1.5604683 

0.568 

-1.683 ± 0.008 

fg2 

— V 2 + 2vj 

1.7086494 

0.160 

-1.929 ± 0.027 

— vi+ 2 v2 

1.8118958 

0.083 

-0.409 ± 0.053 

— I'l + 1^2 + ^7 

1.9853866 

0.135 

-0.672 ± 0.032 

— vi+ 2vj 

2.1588774 

0.094 

-1.110 ± 0.046 

1^1 + 1^2 

2.2731076 

0.294 

-3.055 ± 0.015 

Vl + Vj 

2.4465984 

0.322 

2.872 ± 0.014 

2v2 

2.7233356 

1.244 

-3.020 ± 0.004 

1^2 + ^3 

2.7448194 

0.096 

1.771 ± 0.045 

1^2 + ^7 

2.8968264 

2.492 

2.985 ± 0.002 

V3 + V7 

2.9183102 

0.081 

1.707 ± 0.054 

— Vl + 2 v 2 + 21V3 — vj 

3.0430404 

0.058 

-1.096 ± 0.074 

2v7 

3.0703172 

1.184 

2.628 ± 0.004 

— Vl + 2 v 2 + vj 

3.3470544 

0.430 

-2.237 ± 0.010 

— Vl + V2 + 2 v7 

3.5205452 

0.314 

-2.588 ± 0.014 

2 v 2 + 1^7 

4.2584942 

0.064 

2.436 ± 0.068 

V2 + 2v7 

4.4319850 

0.076 

2.071 ± 0.057 




Figure 8. Phasor plots for for the 7 Dor star KIC 7468196 (top) and for 
the SPB star KIC 10799291 (bottom). It can be seen that the difference be¬ 
tween the light curves of these two stars (Fig. 6) - upward for the 7 Dor 
star and downward for the SPB star - is a consequence of the phases of the 
combination frequencies. 
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Figure 9. A section of the light curve of KIC 10118750 spanning 20 d show¬ 
ing the symmetric, then downward character of the light variations. The 
time is relative to BJD 245 5000. 


6 TWO MORE SPB STARS 

To emphasise the point that the SPB stars are purely g-mode pul¬ 
sators, and to extend the application of combination frequencies 
to understanding these stars, we show two further cases where we 
do not agree with the classification of the star in the literature: 
KIC 10118750, which has a downward light curve and has been 
suggested to be a /3 Cep star, and KIC 5450881, which has a light 
curve that would probably have been called ‘symmetric’ in previ¬ 
ous discussions and which has been previously classified as a rota¬ 
tional or orbital variable. In both cases we view these stars as SPB 
stars with combination frequencies. 


6.1 KIC 10118750 

Fig. 9 shows a section of the light curve of KIC 10118750 where the 
at first symmetric, then downward appearance can be seen. McNa¬ 
mara, Jackiewicz & McKeever (2012) classified this star as a /3 Cep 
star, probably because of the higher frequency range and frequency 
groups compared to some SPB stars. We find that there are only five 
pulsation modes, all with frequencies in the 3.6 — 4.2 d~^ range, 
which could be consistent with a /3 Cep classification. However, the 
combination frequencies, the downward light curve and, especially, 
the late-B effective temperature, T^s = 11400 K, suggest instead 
that this is an SPB star. 

KIC 10118750 has five principal base frequencies that we at¬ 
tribute to g-mode pulsations; most of the rest of the variance is 
in the combination frequencies. Fig. 10 shows the amplitude spec¬ 
trum with frequency groups fgO, fgl and fg2 marked in the top 
panel. There are higher frequency groups of much lower ampli¬ 
tude that can be seen, but are not shown here. The bottom panel 
expands the frequency range around fgl and labels the five base 
frequencies selected. Those five frequencies and their combination 
frequencies with terms up to order 2v were fitted by least-squares to 
the QO-17 data set. We restricted the fit to the five base frequencies 
and only those combination frequencies that had amplitudes greater 
than 30 qmag in the frequency range up to 10 d~^, resulting in 31 
frequencies given in Table 5. 

The middle panel of Fig. 10 shows the stunning removal of all 
peaks in fgO, fgl and fg2 when the 31 frequencies from Table 5 
are fitted to the data. Even the very low amplitude peaks remain¬ 
ing can be modelled with additional combination frequencies®. We 
conclude that KIC 10118750 is an SPB star with a relatively high 
g-mode frequency range for the class. 



frequency d’^ 
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Figure 10. Top panel: An amplitude spectrum for KIC 10118750 out to 
10d~^. There are no p mode pulsations; the higher frequency peaks are 
combination frequencies from five pulsation mode frequencies in fgl. Be¬ 
cause the higher frequency groups are composed of higher order combi¬ 
nation frequencies, we concentrate on the groups fgO, fgl and fg2 here to 
show the result. Middle panel: an amplitude spectrum of the residuals after 
fitting the five pulsation frequencies and 26 combination frequencies with 
amplitudes above 30 p.mag. There are significant peaks remaining with am¬ 
plitudes below our imposed 30 ftmag limit; further combination frequencies 
probably explain those, too. The ordinate scale has been kept the same in 
all panels to give the impact of the variance reduction with such a simple 
explanation. Bottom panel: An expanded amplitude spectrum of the Kepler 
Ql-17 data for KIC 10799291 out to 10d~^ with the pulsation frequen¬ 
cies to 1/5 marked. Note that the lowest unmarked peak in fgl could be 
substituted for i/4 and the same solution would be found. 


6.2 The phases 

The strong downward light curve is described by combination fre¬ 
quencies with relative phases near tt. The minima in this star reach 
abnormal depths, even for stars with asymmetric light curves, as is 
seen to the right of Fig. 9. As such, the phasor plot in Fig. 11 shows 
points that not only lie predominantly on the left, but that also have 
amplitudes of similar magnitude to their base frequencies. 


^ We chose to keep the number of fitted frequencies low for the impact 
this has. It is not our purpose here to generate thousands of combination 
frequencies to explain all of the tiny variance remaining. 
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Table 5. A least-squares fit of the five base frequencies of KIC 10118750 
and the combination frequencies with terms up to order 2z/, with frequen¬ 
cies up to 10d~^, and with amplitudes restricted to greater than 30 p-mag. 
There are 31 identified frequencies, including the five base frequencies. The 
zero point of the time scale is BJD 245 5694.25. 


labels 

frequency 

d-i 

amplitude 

mmag 

± 0.002 

phase 

radians 

fgO 

— 2 u 2 + I'z 

0.0069965 

0.038 

2.279 ± 0.059 

— 2 u 2 -\- 2 i /'4 — lx's 

0.0077650 

0.066 

1.692 ± 0.033 

2 i^l — 21^3 — i/'4 -|- 1^5 

0.0136230 

0.032 

-1.947 ± 0.070 

— I'l + 1^2 

0.0473818 

0.309 

-2.138 ± 0.007 

-i/l + — t'a + 2 i^4 - 

0.0481503 

0.026 

-2.966 ± 0.087 

-^2 + *^3 

0.0543783 

0.046 

0.403 ± 0.048 

V 2 — 2 iv 3 — M 4 + z /5 

0.0610048 

0.040 

-1.670 ± 0.056 

— 2 i^l + 21^2 

0.0947636 

0.038 

2.440 ± 0.059 

— -1-1^3 

0.1017601 

0.113 

-1.050 ± 0.020 

— - 1 - z^4 

0.3196718 

0.155 

-1.740 ± 0.014 

— 2 v\ -h 1^2 ^4 

0.3670536 

0.052 

2.598 ± 0.042 

fgi 

21^1 — 1^4 

3.3105639 

0.655 

-1.634 ± 0.003 

2 ui — 1^3 + L /4 — 1^5 

3.3113324 

0.030 

-2.996 ± 0.075 

I'l + 1^2 — 1^4 

3.3579457 

0.058 

-0.470 ± 0.038 

— 1^3 + ^4 — ^5 

3.3587142 

0.034 

-1.589 ± 0.066 

1^2 + 1^4 - 1^5 

3.4604743 

0.066 

0.647 ± 0.034 


3.6302357 

1.645 

-1.213 ± 0.001 

^2 

3.6776175 

1.581 

-2.979 ± 0.001 

^3 

3.7319958 

0.570 

-1.726 ± 0.004 

1^1 — 1^2 + 1^4 

3.9025257 

0.049 

-0.597 ± 0.045 

1/4 

3.9499075 

0.796 

-2.697 ± 0.003 

^5 

4.1670507 

0.601 

-2.628 ± 0.004 

fg 2 

I'l + 1^2 + ^4 — ^5 

7.0907100 

0.141 

2.409 ± 0.016 

ui + 2^2 — U 3 

7.2534749 

0.051 

-1.693 ± 0.043 

2 u\ 

7.2604714 

1.151 

-1.624 ± 0.002 


7.3078532 

0.131 

2.503 ± 0.017 

2^2 

7.3552350 

0.401 

0.513 ± 0.006 

I'l + 1^3 

7.3622315 

0.196 

-2.437 ± 0.011 

+ V4 

7.5801432 

0.894 

2.767 ± 0.003 

U 2 + i^4 

7.6275250 

0.047 

0.670 ± 0.047 

t'l + 1^5 

7.7972864 

0.036 

2.864 ± 0.062 



Figure 11. The phasor plot for the SPB star KIC 10118750. 
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Figure 12. A section of the light curve of KIC 5450881 spanning 10 d. The 
time is relative to BJD 245 5000. Almost all of the variation in the light 
curve is explained by two base frequencies and six combination frequencies. 
A rotational hypothesis with spots is not viable. This is pulsation. 
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6.3 KIC 5450881 - the simplest case of a frequency group 
star 

McNamara, Jackiewicz & McKeever (2012) listed KIC 5450881 as 
either a rotational or binary variable. Fig. 12 shows a 10-d section 
of the light curve, which can be described with two base frequen¬ 
cies and six combination frequencies. This can be seen in Fig. 13 
where a dominant triplet and combination frequencies are evident. 
Because the triplet is equally spaced, we only need two base fre¬ 
quencies; we have arbitrarily chosen the two labelled. We then find 
only six combination frequencies, one of which is the other mem¬ 
ber of the triplet. With the two base frequencies and six combina¬ 
tion frequencies, almost all of the variance is strikingly removed, as 
can be seen in the lower panel of Fig. 13. At higher amplitude res¬ 
olution additional significant peaks can be seen in this frequency 
range and at higher frequency; these can be modelled in a more 
detailed analysis in the future. 

This pattern of two base frequencies, which we propose are 
from two g modes, and six combination frequencies gives an ele¬ 
gant explanation of the light curve that is preferable to a rotational 
spot model or a binary star model. Spots and binary variations do 
not give rise to combination frequencies. This star is thus the sim¬ 
plest case of frequency group star, with only two base frequencies. 
The other stars presented in this paper with more base frequencies 
are more complex examples of the same physics, hence this star 
illustrates the principles most clearly. 

The equally spaced frequency triplet could have one of sev¬ 
eral origins. It could be interpreted as a rotational triplet, hence 
probably arising from a g-mode dipole. Alternatively, it could arise 
from pure geometric amplitude modulation, such as is seen in roAp 
stars, or pure frequency modulation, as could arise if there were a 
very massive companion (Shibahashi & Kurtz 2012). Both of the 
latter possibilities can be tested. We chose a time, to, that set the 
phases equal for the two sidelobes to the central peak. With two si¬ 
nusoidal variations this is always possible. We then expect for pure 
amplitude modulation that the central frequency should have the 
same phase at to, and for pure frequency modulation the central 
frequency should have a phase 7t/2 different to the sidelobe phase. 

Neither of these conditions is obtained in this case - the 
phase difference between the central frequency and the sidelobes 
is 2.2 rad - hence we conclude that the triplet arises neither from 
pure amplitude modulation, nor from pure frequency modulation. It 
can then be represented with a rotational triplet, or simply two base 
g-mode frequencies with combination frequencies. If the triplet is 
from rotationally split dipole modes of of order m = —1,0, 4-1, 
then the the inverse splitting of the triplet is l/(iZ 2 — t^i) = 
2.4466 d, and an estimate of the Ledoux constant leads to a ro¬ 
tational period for the star somewhat less than twice that, i.e., 
~ 4 — 5 d. A spectroscopic measurement of v sin i could constrain 



Figure 13. Top panel: An amplitude spectrum for KIC 5450881 out to 
10 showing a triplet and combination frequencies. The two marked 
peaks were chosen as the base frequencies to generate the combination fre¬ 
quencies; the other member of the triplet would serve equally well. Bottom 
panel: the amplitude spectrum of the residuals after pre-whitening by the 
two base frequencies and six combination frequencies with terms up to or¬ 
der 2u. 


Table 6. A least-squares fit of the two base frequencies of KIC 5450881 
and six combination frequencies with terms up to order 2v, where one of 
the combinations is the other member of the dominant triplet. The zero point 
of the time scale is BJD 245 5694.25 


labels 

frequency 

d-i 

amplitude 

mmag 

±0.001 

phase 

radians 

— 1^1 + 1^2 

0.4087346 

0.110 

0.9693 ±0.0107 

— 2^1 + 21/2 

0.8174692 

0.007 

1.4095 ±0.1584 


3.2699899 

0.768 

1.4283 ±0.0015 

^2 

3.6787245 

1.117 

2.5033 ±0.0011 

-Ui 2^2 

4.0874591 

0.779 

-0.8178 ±0.0015 

2vi 

6.5399798 

0.013 

0.0892 ±0.0921 

yi + V2 

6.9487144 

0.103 

-0.4452 ±0.0115 

2V2 

7.3574490 

0.886 

-2.6611 ± 0.0013 


or refute that possibility. Breger & Kolenberg (2006) examined sev¬ 
eral cases of nearly equally split triplets in a variety of pulsating 
stars and suggested that in those cases the triplets were better un¬ 
derstood in terms of two base mode frequencies and combination 
frequencies. That may be the case, also, for KIC 5450881. 


6.4 The phases 

The light curve of KIC 5450881 is not strongly asymmetric, hence 
the combination frequencies describing it are few in number and 
are not heavily biased to one side of the phasor diagram shown in 
Fig. 14. Note that a large amplitude hut a relative phase close to 
7 t/2 or 37t/2 does not imply an asymmetric light curve; that is, a 
symmetric light curve could have many such peaks. 


7 Be STARS 

Several types of stars are given the spectral classification Be (B 
emission line) stars. There are pre-main-sequence Herbig Be stars 
(usually grouped with the Ae stars), there are binary mass transfer 
systems where circumstellar emission occurs, and there is a large 
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Figure 14. The phasor plot for the SPB star KIC 5450881. 


class of stars that are rapidly rotating and have occasional outbursts 
where material is launched from the stellar surface into a disk, thus 
producing emission lines. Members of the latter group have been 
shown to be pulsating stars (Rivinius, Baade & Stefl 2003) and a 
model has been proposed by Dylan Kee et al. (2014) of how pulsa¬ 
tion, coupled with rapid, but sub-critical rotation, can launch mate¬ 
rial into orbit. 

During the outbursts of pulsational Be stars the pulsation am¬ 
plitude increases. Some B stars observed by Kepler show pulsa¬ 
tional outbursts at long time intervals, and we consider those to be 
pulsational Be stars, even though we do not have spectroscopic ob¬ 
servations to show the line emission that should accompany these 
pulsation outbursts. Future spectroscopic observations can test this. 

We show here one example for which the light curve and am¬ 
plitude spectrum - during and outside of pulsational outburst - are 
strikingly similar to the frequency groups of SPB stars and also 
the 7 Dor stars that we have presented in this paper. This similarity 
gives strong support to the hypothesis that the pulsational Be stars 
eject material when several, or many, g modes come into phase with 
each other. We suppose that those modes are prograde sectoral g 
modes where the pulsation velocity is primarily horizontal and in 
the direction of rotation, so that it adds to the subcritical rotation 
velocity to reach critical velocity and launch a circumstellar disk. 

In a review of pulsation in B stars and the prospects of as- 
teroseismology, Aerts (2006) concluded: ‘It seems that pulsating 
Be stars are complicated analogues of the SPB stars ..We con¬ 
cur: The single pulsating Be stars are rapidly rotating SPB stars. 
No physics other than pulsation need be conjectured to understand 
these stars. 


7.1 An example of a pulsating Be star: KIC 11971405 

KIC 11971405 (HD 186567) is a relatively bright (V = 9.2) B8 V 
star that shows g-mode pulsations typical of SPB stars. McNamara, 
Jackiewicz & McKeever (2012) listed it as a frequency group star, 
while Pigulski et al. (2009) describe the star as aperiodic from 
ASAS ground-based data. We show here that it is an SPB star with 
frequency groups, and with occasional pulsation ‘outbursts’ typical 
of pulsating Be stars. We hypothesise that it is a Be star, and that 
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Figure 15. Top panel: The full Kepler QO-17 light cuive of KIC 11971405. 
At this scale the pulsations are not visible, but the pulsation outbursts later in 
the mission can be seen. Second and third panels: 20-d sections of the light 
curve where two of the outbursts occur. The variations are typical symmet¬ 
ric or downward SPB pulsations. Fourth panel: a 20-d section outside of 
outburst at the same scale as the previous light curves. Bottom panel: the 
same as the fourth panel, but at expanded scale for visibility. 


the pulsation outbursts are the result of modes with very closely 
spaced frequencies that only come into phase at long intervals to 
build up to an outburst. 

Fig. 15 shows light curves for KIC 11971405. The top panel 
shows the full Kepler QO-17 data set. At this scale it is not possible 
to see the individual pulsations, but the pulsation outbursts later 
in the mission are evident. In the second and third panels, 20-d 
sections of the light curve show two of those outbursts in detail. 
The fourth panel shows a typical 20-d section of the star during 
quiescence at the same scale as the outburst light curves, and the 
bottom panel expands the ordinate scale for a more detailed look 
at the quiescent pulsational variations. KIC 11971405 looks like a 
typical SPB star, except for the occasional outbursts. 

Fig. 16 shows amplitude spectra for the entire QO-17 data set 
(top) and for the independent subsets of the data QO-9 and QIO- 
17 (middle and bottom). There are three main frequency groups, 
with further groups at much lower amplitude at higher frequency 
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Figure 16. Top panel: An amplitude spectrum for KIC 11971405 for the en¬ 
tire QO-17 data set out to 6d~^. There are some higher frequency groups 
at much lower amplitude that are not shown. Frequency groups fgO, fgl and 
fg2 are marked. Although the highest amplitude peaks are in fg2, those are 
combinations of frequencies in fgl. Middle and bottom panels: The same 
as the top panel, but for independent data subsets QO-9 and QlO-17. These 
show that the amplitudes of some of the peaks have changed, but the fre¬ 
quencies and the frequency patterns remain fixed. 
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Figure 17. Top panel: An amplitude spectrum for KIC 11971405 for the 
QO-9 data set with some peaks in fgl labelled. We use these as base fre¬ 
quencies to calculate combination frequencies and show that the high am¬ 
plitude peaks in fg2 and fgO are simple combinations of frequencies in fgl. 
Middle panel: an amplitude spectrum for QO-9. Bottom panel: An ampli¬ 
tude spectrum for QO-9 after pre-whitening the 15 frequencies selected in 
fgl and four combination frequencies that remove most of the variance in 
fgO and fg2. 


that are not displayed at this scale. This figure shows that the am¬ 
plitudes of some of the peaks change between the independent data 
subsets, but the frequencies and their patterns remain the same. We 
hypothesise that this is caused by some pulsation mode frequencies 
that are so close in frequency that they are not resolved in the 4- 
yr time span of the Kepler data set. Those beat against each other 
slowly - hence the long time span between pulsation outbursts - 
and also come into phase at times with other pulsation modes. Out¬ 
bursts arise when many modes come into phase to give the larger 
amplitudes seen in the light curves in the second and third panels 
of Fig. 15. 

A closer look at the amplitude spectrum for QO-9, where there 
were no large outbursts, shows that the amplitude spectrum is dom¬ 
inated by combination frequencies. In particular, the highest am¬ 
plitude peaks in fg 2 are all simple combinations of frequencies in 
fgl. In Section2 we showed how it is possible for the combination 
frequencies to have higher amplitudes than the base frequencies. 

Fig. 17 shows the amplitude spectrum for the QO-9 data in 
the range of fgl with a series of peaks identified. We selected from 
these our base frequencies for calculating the combination frequen¬ 
cies. Because there are so many frequencies in this group, we chose 
not to generate the thousands of combination terms that possibly 
arise among them, but instead selected just a few simple combina¬ 
tions to illustrate that the peaks in fg2, in particular, but also fgO, 
are combination frequencies. The frequency of the peak in fgO is 
equal to the difference of the frequencies of the highest amplitude 
peaks in fg 2 , but those themselves are simple sum combination fre¬ 
quencies of base frequencies in fgl. 

We have thus shown that the SPB star KIC 11971405, which 
shows pulsational outbursts, has a light curve and amplitude spec¬ 


trum similar to those of other SPB stars and 7 Dor stars that show 
frequency groups. The principal difference is that the amplitudes 
of some of the peaks in the amplitude spectrum of KIC 11971405 
are variable over the 4-yr time span of the data set. We hypothesise 
that this is caused by pulsation modes with frequencies too closely 
spaced to be resolved in 4 yr. This example and others in the Kepler 
data, support the hypothesis that the pulsational Be stars are rapidly 
rotating SPB stars. 


8 DISCUSSION AND CONCLUSIONS 

The frequency spectra of the g-mode pulsators of the main- 
sequence, the 7 Dor stars and the SPB stars, show a wide variety 
of complexity. Some stars show series of consecutive radial over¬ 
tone g mode frequencies, and it is those that are the most interest 
asteroseismically, since they allow inference about the core condi¬ 
tions in the stars (e.g.. Van Reeth et al. 2015; Saio et al. 2015; Kurtz 
et al. 2014; Bedding et al. 2014; Bouabid et al. 2013). Other stars 
show frequency groups that have not been understood, or have been 
only partially understood, until now. We have shown in this paper 
that the frequency groups found in the 7 Dor stars and SPB stars 
are often dominated by combination frequencies of only a few base 
frequencies, and we interpret those base frequencies as arising from 
a small number of g modes. 

We have shown that the shapes of the light curves, which have 
previously been classified by visual inspection into ‘symmetric’ 
and ‘asymmetric’ are a consequence of the same physics. It is the 
phases of the combination frequencies that describe the visual ap- 
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Table 7. A least-squares fit of the 15 frequencies in fgl of KIC 11971405 
plus four combination frequencies from four base frequencies in the fgl 
set. The three highest peaks in fg2 and the peak in fgO are thus shown to 
be combination frequencies. The residual amplitude at those frequencies 
is a consequence of the change in amplitude, probably because unresolved 
frequencies close to the base frequencies. The zero point of the time scale 
is BJD 245 5694.25 


labels 

frequency 

d-i 

amplitude 

mmag 

±0.004 

phase 

radians 

fgO 

1^4 + 1^15 -I'i- I'M 

0.276260 

1.858 

0.809 ± 0.002 

fgl 

1^1 

1.779808 

0.236 

-2.532 ±0.017 

^2 

1.785546 

0.318 

0.200 ±0.013 


1.801339 

0.207 

-0.644 ± 0.019 

1/4 

1.822719 

0.859 

1.824 ± 0.005 

1^5 

1.838304 

0.736 

0.303 ± 0.005 

UQ 

1.863744 

0.383 

0.451 ± 0.010 

1/7 

1.871133 

0.163 

-2.767 ± 0.025 


1.879095 

0.362 

0.541 ± 0.011 


1.887687 

0.430 

2.551 ± 0.009 

^10 

1.895767 

0.850 

3.116 ±0.005 

^11 

1.915925 

0.309 

-0.289 ±0.013 

^12 

1.948357 

0.227 

2.728 ±0.018 


2.094494 

0.287 

2.116 ±0.014 

VXA 

2.171925 

1.166 

0.071 ± 0.003 


2.187612 

1.826 

-1.173 ±0.002 

fg2 

!Z5 + 1^10 

3.734071 

1.571 

-2.086 ±0.003 

I'S, + Viz 

3.973589 

1.412 

-0.441 ± 0.003 

V4 + VIS, 

4.010331 

3.867 

1.032 ±0.001 


pearance of the light curves. We have extended this to show that 
‘downward’ light curves, which some investigators have conjec¬ 
tured are caused partially or completely by stellar spots, are also 
explained by a small number of g mode pulsations and their com¬ 
bination frequencies. We further showed that SPB stars with pulsa- 
tional outbursts show frequency groups and combination frequen¬ 
cies that are the same in character to other SPB stars, supporting 
to the conclusion that pulsational Be stars are rapidly rotating SPB 
stars. 

The only physics necessary to understand the light curves of 
the 7 Dor stars, the SPB stars and the pulsating Be stars is g-mode 
pulsation. For many years, there has been conjecture that spots are 
a source of the light curve shapes in these stars. A particular driver 
for that conjecture is the shape of the light curves of the stars that 
show ‘downward’ light curves. We have shown that pulsation alone 
accounts for the light curve shapes. 

For the 7 Dor and SPB stars presented in this paper, the base 
frequencies and the combination frequencies are unchanged to high 
precision over the 4-yr time span of the data. No spot model can ac¬ 
count for this. We know what stable spots look like on upper main- 
sequence stars from the chemically peculiar magnetic Bp, Ap and 
Fp stars. Those show rotational light curves that are non-sinusoidal, 
hence can have many harmonics of the rotational frequency. But 
there are no combination frequencies generated. In the frequency 
groups of the stars studied in this paper, the combination frequen¬ 
cies dominate much of the amplitude spectrum and harmonics of 
the base frequencies have low amplitudes; this is not the signature 
of spots. 


For asteroseismic studies of SPB stars and 7 Dor stars, combi¬ 
nation frequencies have been considered a nuisance to be detected 
and discarded in the acquisition of pulsation mode frequencies for 
asteroseismology. We now see that in some cases they can com¬ 
prise the major components of the amplitude spectra. Indeed, for 
asteroseismology they may still be modelled and removed. But it 
may be possible to use the combination frequencies for asteroseis¬ 
mic study themselves, as has been done for white dwarf stars for 
many years and as is being attempted for p mode pulsations in 5 Set 
stars. The large mode energies implied by the visibility of surface 
g modes promise the discovery of interesting new physical infor¬ 
mation about the cores of these main-sequence stars. 
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